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Within the Free- Volume Asakura-Oosawa-Vrij (FVAO) theory of colloid-polymer mixtures, we 
show that unphysical gas-liquid binodals predicted in the regime of small attraction range (i.e. 
polymer size) are caused in part by the use of the Carnahan-Starling (CS) hard sphere (HS) reference 
free energy. Replacement of the CS expression with an alternative dramatically affects predicted 
phase behaviour and, for polydisperse colloid, the resultant fractionation predictions. Although 
short-range attractions render FVAO, as a perturbative HS-based theory, less accurate anyway, we 
argue that the particular effects of CS in this regime are an important consideration - usually ignored 
- in the evaluation of such theories. We refer to a variety of literature exhibiting similarly inaccurate 
gas-liquid binodals, and suggest CS's status as the de facto choice of hard sphere reference should 
be carefully considered where short-range attractions are present. 

PACS numbers: 64.10. +h,61. 20. Gy,82.70.-y,82.70.Dd,64.75.-g 



I. INTRODUCTION 

Long ago, Widom [l| explained why excluded-volume, 
or equivalently entropic, effects should dominate the be- 
haviour of liquids around the triple point. The long-range 
inter-particle attractions (e.g. as modelled by a van der 
Waals potential) of > 10 nearest neighbours sum approxi- 
mately to a constant. To a first approximation, therefore, 
the inter-particle attraction simply lowers the internal en- 
ergy component of the Hclmholtz free energy by a fixed 
amount. Minimising the free energy is then equivalent 
to maximising the entropy. This qualitative argument 
underpins the success of the use of repulsive particles in 
general , and hard spheres in particular || , as reference 
systems in perturbation theories of the liquid state. 

In the hard sphere (HS) model, particles interact via 
a step-like repulsive potential which is infinite at contact 
and zero everywhere else: 
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where d = 2a is the HS diameter and r the distance 
between particle centres. This model has a long his- 
tory in the study of fundamental statistical mechanical 
problems of the condensed state. Thus, e.g., there was 
the celebrated question of whether a collection of hard 
spheres would crystallise at high enough number density 
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[J], which was eventually settled in the affirmative first 
using computer simulations @, and then more recently 
using experiments on colloids in which the particles in- 
teract as nearly-perfect hard spheres 0. These latter 
experiments took hard spheres into the laboratory. They 
have now been used to study everything from the kinet- 
ics of crystal nucleation Q to the glass transition Q. In 
such experiments, very direct confrontation with analytic 
calculations and simulations is possible - there is in prin- 
ciple only a single, experimentally accessible, parameter - 
volume fraction Q , which for N hard spheres of diameter 
d filling a volume V is <j> = nNd 3 /6V. 

Inter-particle attraction can be induced in colloids 
by adding non- adsorbing polymers: exclusion of poly- 
mer molecules from the region of space between the 
surfaces of two nearby particles results in a net os- 
motic force pushing the particles together. The range 
of this 'depletion' attraction is controlled by the size of 
the polymer, e.g. as measured by its radius of gyration 
r g . The strength of the depletion attraction is propor- 
tional to the polymer's chemical potential, or, equiva- 
lently, the number density tir of polymers in a reser- 
voir in osmotic equilibrium with the colloid-polymer mix- 
ture. A well-characterised colloid-polymer (CP) mixture 
of this kind, hard-sphere-like stcrically-stabilised parti- 
cles of poly-methylmcthacrylatc (radius a) and the non- 
adsorbing polymer polystyrene, has been used to study 
a range of basic condensed matter questions (reviewed in 
Ref. |l0|). In particular, experiments using this system 
have demonstrated that a gas-liquid coexistence region 
only exists in the equilibrium phase diagram when the 
depletion attraction is of long enough range: r g /a > 0.24. 
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A simple model of the depletion attraction, the 
Asakura-Oosawa-Vrij (AO) model [ill, EH, treats the 
colloids as hard spheres, of radius a; the polymers are 
(smaller) spheres, of radius S « r g that interpenetrate 
one another, but have a hard sphere interaction with the 
colloids. The AO model has been incorporated into a 
free- volume theory for calculating the phase behaviour of 
CP mixtures [Sj. This Free- Volume AO (FVAO) theory 
requires as input reference free energies of hard spheres 
in the fluid and crystalline states. For the latter, a pa- 
rameterised expression from simulations [bij is available. 
For the HS fluid, a simple closed form expression given 
by Carnahan and Starling (CS) [l5| is the default choice 

[31 Ml- 

The CS free energy is obtained from integrating the CS 
equation of state (EOS), which is usually given in terms 
of the compressibility factor, Z = PV/mkBT relating the 
pressure P, volume V and temperature T of m moles of 
hard spheres: 
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This is known to represent simulation data accurately, 
at least at low to moderate 0, but is self evidently un- 
physical at high <fi, since it does not diverge until 0=1, 
whereas a monodispcrse hard sphere fluid may only be 
packed to random close packing, 0rcp ~ 0.64, and 
even the ordered (crystalline) close packing value is only 
<?Wc = 7r/3\/2 « 0.74. While there is by no means con- 
sensus as to the "correct" behaviour of the hard sphere 
equation of state at high 4> EH (especially given that, 
very close to 0rct>, jamming and path-dependent effects 
come into play [201), it is clear that CS's high <p behaviour 
cannot be accurate. This unphysical feature of the CS 
EOS matters little in many contexts. However, we argue 
in this note that it contributes to unphysical predictions 
of phase behaviour in the FVAO theory of CP mixtures 
in the regime of small size ratios, q = 5/a -C 1. On re- 
placing the CS EOS with one that is more realistic at 
high </>, we find a marked change in the predicted phase 
behaviour at small q. This leads also to important conse- 
quences when the colloid in a CP mixture is polydisperse: 
the fractionation upon phase separation predicted by the 
two approaches differ by an order of magnitude. 

Unphysical gas-liquid binodal predictions in FVAO 
and those in similar theories PU are usually attributed 
to the increasing unsuitability of HS-based perturbative 
theories as attractions become shorter-ranged: correla- 
tions become less HS-likc, and the HS reference itself be- 
comes significantly mctastable. We stress from the out- 
set that this remains the case no matter what choice of 
HS reference is used - modifying FVAO by replacing CS 
still results in poor comparisons with available simula- 
tion data, and recent work improving on e.g. the approx- 
imation of non-interacting polymers [2lj remains neces- 
sary and well-motivated. However, we argue that the 
tendency of shorter-range attractions to bring particles 



closer together necessarily renders the high <f> behaviour 
of a HS reference more important, in a variety of theo- 
ries, so that testing specifically for dependence on CS is 
an important step in evaluating their predictions. 



II. THE FVAO THEORY 



A. Free energy and hard sphere compressibility 
factor 



In treatments of CP mixture phase behaviour prior to 
the FVAO theory, the AO depletion potential (which, for 
q < 0.154, maps exactly to an effective pair potential |18j ) 
was used as a perturbation to the HS reference. In FVAO, 
the free energy contribution of the (ideal) polymers is 
obtained directly, explicitly allowing for partitioning of 
polymers between the phases. Within FVAO [l3[, the 
Hclmholtz free energy per unit volume of the mixture is 
found to be the sum of two terms (with fceT = 1): 
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where Zjjs(0) is the HS compressibility factor, a the col- 
loid radius, n the polymer number density and tir the 
polymer number density in a reservoir in osmotic equi- 
librium with the CP mixture (i.e. coupled to the system 
via a membrane through which only polymer and solvent 
may pass). Note that the apparent separation in colloid 
and polymer terms in this expression is deceptive, since 
n-n is related to the concentration of polymers in the vol- 
ume left accessible to them by the colloids: tir = n/a 
where a(4>) is the free volume fraction left by the col- 
loids. We follow the original FVAO prescription [l]| that 
a takes the same value as at zero polymer fugacity and is 
given by a scaled-particle (or equivalcntly, Pcrcus-Yevick) 
[22| result. We note that this assumption - which renders 
FVAO, too, an essentially perturbative theory - becomes 
less appropriate as attraction range decreases [l8], but it 
is sufficient for our purpose of comparing an approximate 
treatment involving CS to another choice of reference. 

To apply Eq. [31 fluid and crystal branches of Zhs are 
required, with integration constants chosen to reproduce 
the known fluid-crystal coexistence gap in the HS phase 
diagram. Since we only deal with gas-liquid coexistence 
in this work, integration constants are not required. 

In the original implementation of the FVAO theory and 
most subsequent applications, Zjjs = -^cs, Eq. [5J is used 
as the fluid compressibility factor in Eq. [3J We compare 
Zcs with an alternative approximate form, the 'WKV3' 
expression originally given in Ref. (23j . It is given by: 
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where {Ih 



denote a set of seven coefficients, and £ = 
. Note that this expression does not diverge 
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FIG. 1. The difference, Zwkv3 — Zcs> between the two HS 
fluid compressibilities considered, as a function of volume frac- 
tion. Inset: behaviour of Zwkv3 (dashed) and Zcs (full line). 



at </>rcp, as might be expected for a HS fluid, but rather 
at c i osc . A recent review [l9j recommends Z\vkv3 as 
the best expression for the HS fluid at high <j>, but the 
difficulty of simulating in this region to determine what 
constitutes "best" qualifies such a recommendation. In 
the present context, the important feature of WKV3 is 
that it closely matches CS at low to moderate <f>, but 
differs on approaching 0rcp and above. This difference 
should become more important as attraction range de- 
creases, so WKV3 serves as a good choice to illustrate 
the sensitivity of predictions to the choice of HS refer- 
ence. Fig. Q] shows that Zcs and Z WKV 3 are nearly iden- 
tical for < 0.55, but Zwkv3 overtakes Zcs beyond that 
point, with the difference becoming increasingly marked 
as <j) — > 0R,cp • We now turn to explore the consequences 
of this difference for phase behaviour at small q. 



B. Phase behaviour 

A full implementation of the FVAO theory encompass- 
ing fluid and crystal phases shows that gas-liquid coex- 
istence with monodispcrse colloid is thcrmodynamically 
stable only when q > 0.3. At low q, therefore, the un- 
physical divergence of Zcs as cj) — > 1 has little effect on 
the equilibrium phase diagram, since the density of the 
crystal phase is determined largely by the behaviour of 



the Hall expression [14J near 



0.74. However, the 



metastable gas-liquid binodal remains important, since it 
plays a key role in phase transition kinetics [24l - l33j . e.g. 
enhancing crystal nucleation [25L I30T ] , splitting advanc- 
ing crystal interfaces [3l|, |32T|. o r causing nonequilibrium 
behaviour such as gelation [33[. Here we therefore ne- 
glect the crystal branch of the free energy and consider 
the gas-liquid binodals predicted by Eq. [3] with Zcs and 

Z\VKV3- 



FIG. 2. FVAO gas-liquid binodals calculated with two choices 
of Zhs, Eqs.[2]and[4] The filled symbols show Zcs, the open 
ones Zwkvs for q = S/a = 0.4(o, •), 0.15(0, ♦), 0.1(D, ■). 
Dashed line: (f> — 0.64 « </>rcp- For comparison, the solid 
line shows the simulated binodal for q — 0.15 from Ref. [30l ] . 
and the dotted line that for q = 0.1 from Ref. [3]. At q = 0.1, 
recent simulations [34|] find a critical colloid packing frac- 
tion <j> = 0.249. However, the position of the critical point 
(where mean-field theories are unsuitable) is not the focus of 
the present work, rather the position of the liquid binodal at 
higher nu- 



By equating relevant derivatives of Eq. [31 which 
give the pressures and chemical potentials in coexisting 
phases, we calculate the gas-liquid coexistence bound- 
aries (or binodals) at different q [l3[ with Zrs given by 
Eqs.[5]and[3J Fig. [5] Note that the vertical axis here gives 
the polymer reservoir concentration, which is directly re- 
lated to the polymer chemical potential, so that tic lines 
in this representation are horizontal. 

At the moderate size ratio of q = 5/a = 0.4, the two 
choices of Zhs give similar results. However, when q is 
reduced, utilising Zcs results in the liquid side of the bin- 
odals exceeding </>rcp > i- e - t ne maximum possible volume 
fraction of a HS fluid, which is unphysical. At q = 0.1, 
even the maximum crystalline packing, </> c iose ~ 0.74, is 
exceeded. This is because when Zrs is used to calculate 
coexistence, the liquid does not pay a high enough free 
energy penalty for approaching and then exceeding un- 
physical concentrations. In contrast, in the same range of 
q, utilising Zwkv3 results in binodals whose liquid sides 
stay below ^rcp- 

Fig. [2] includes simulated gas-liquid binodals at low- 
q from Refs. [TH, which still compare poorly with 
FVAO when WKV3 is used (note also that the q = 0.1 
critical colloid packing fraction </> = 0.249 found in re- 
cent simulations [3~4| is approximately a factor of 2 lower 
than the apparent location using FVAO and WKV3). As 
discussed in the Introduction and extensively in Ref. [l8| 
FVAO, as a perturbative theory assuming zero polymer 
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fugacity, remains inaccurate at \ow-q for reasons addi- 
tional to the use of CS. 

Given this, it is important to emphasise that our aim 
is to illustrate the sensitivity to choice of HS reference 
brought about by short-range attractions, and not to im- 
ply that all inaccuracies in this regime are due solely 
to the use of CS. FVAO simply serves as a well-known 
archetype of a wide class of theories which build upon 
the hard sphere reference. 

We note that unphysically high liquid densities have 
also been found in the gas-liquid binodals predicted by 
more sophisticated models of CP mixtures at low q (e.g. 
2l| ) as well as liquid-state theories of other systems with 
short-range attraction (e.g. [i~7l Il8l. l2ll. 135)439 ] ) . In each 
case cited, the CS reference free energy is used exclu- 
sively and without question. In addition to other dangers 
that can affect any HS reference in this regime [l8j ]. we 
suggest that CS is a specific, important factor whose in- 
fluence on predicted gas-liquid binodals should be tested 
in the course of evaluating predictions at short attrac- 
tion range. Motivating this is the generic observation 
that shorter-range attractions, however they are incor- 
porated into a theory, should act to compress the liquid 
phase, favouring the minimisation of interaction energy 
- broadly speaking, the role of the hard sphere reference 
is to define the penalty the system pays for this compres- 
sion which, for CS, becomes spuriously low at higher </>. 
Thus, as attraction range is decreased within perturba- 
tive (and perhaps even non-perturbative [l3|) theories, it 
becomes more important to test for dependence on CS 
in order to distinguish its effects from those intrinsic to 
the theory in question. 

Concluding this section, we note that for the purpose of 
precisely predicting a critical point, it appears sufficient 
and preferable to use a criterion based on the second 
virial coefficient B2 [H], H(| , which gives good agreement 
down to \ow-q and avoids the concerns outlined above. 
However, the behaviour of binodals away from the crit- 
ical point cannot be discerned by these means, and this 
behaviour is important e.g. for determining the relation- 
ship between binodals and features such as gelation or 
vitrification lines fl7j . predicting crystal growth kinetics 
[24l [29j , andpredicting fractionation between separating 
phases 0, S2|. In respect of the latter, we now turn to 
illustrate the effects of the choice of Zhs 011 fractiona- 
tion predictions in the case of polydisperse colloid, where 
the metastable gas-liquid binodal is particularly impor- 
tant since crystallisation is suppressed in experimental 
realisations. 



III. POLYDISPERSITY AND FRACTIONATION 

When inter-particle attraction causes a gas-liquid type 
transition in a colloidal suspension in which the particle 
size distribution has finite polydispersity (defined as the 
standard deviation divided by the mean of the size distri- 
bution), the two daughter phases not only differ in their 



particle number densities, but also in their size distribu- 
tions. 

We distinguish between two well-known formulations: 
a perturbative approach [43j , and a moment free energy 
approach [44Tj46j . In the former, only a free energy for the 
monodisperse (single-sized) reference system is required, 
whereas the latter requires a polydisperse model free en- 
ergy, such as the Boublik-Mansoori-Carnahan-Starling- 
Leland (BMCSL) one [47|, which generalises CS. In the 
context of CP mixtures, both can be applied to the FVAO 
theory [H, [45J . Since the perturbative approach allows 
substitution of the (monodisperse) ^wkv3 into the ex- 
isting theory, we examine fractionation in this context. 
However, as the BMCSL free energy inherits the 0=1 
divergence of the parent CS expression, we expect the 
findings to be also of relevance to the moment free en- 
ergy approach. Indeed, our findings are likely generic 
to any treatment of fractionation in which some form of 
approximate HS-like free energy is used. 

In the perturbative treatment [43| . the fractionation, 
or difference in mean particle diameter ((e)) between co- 
existing gas (g) and liquid (I) phases, is given by: 



im = -[A/p]y P , 



(5) 



where A(p) = p[dp cx (e)/de] describes, in a monodisperse 
reference system, the variation in excess chemical poten- 
tial /x ex as a function of fractional size deviation e (pos- 
itive or negative) from the overall mean diameter. It 
is the rate of change of excess chemical potential with 
particle size. The variance of the parent distribution is 
(Tp. Thus, —[A/p] g controls the 'fractionation strength' 
- it gives the fractional difference in mean diameter be- 
tween phases, normalised by the overall variance. For 
a CP mixture, plotting —[A/p] l g as a function of poly- 
mer reservoir volume fraction is therefore a convenient 
way of seeing how much fractionation results at various 
strengths of the depletion attraction (see [43[ for details) . 
We show such a plot at different size ratios, q = S/a in 
Fig. [3] calculated within the perturbative framework [43| 
using Z C s and Z WKY3 . 

As in the case of the binodals, Fig. there is little 
difference at q = 0.4; but the two compressibility fac- 
tors give progressively more divergent predictions for the 
strength of fractionation as q decreases. In the case of the 
binodal, the unphysical nature of the predictions made 
using Zcs is immediately obvious on physical grounds: a 
liquid cannot pack at <fi > </>rcp ■ No such intuition guides 
our evaluation of the predictions presented in Fig. [3] be- 
cause to date, there are relatively few detailed experimen- 
tal, theoretical or simulation studies of fractionation. 

Some precise results are emerging for pure HS sys - 
tems (e.g. [48U50l | ) . For CP mixtures, experiments (5lT - 
[53[ have tested and confirmed the scaling of fractiona- 
tion strength with parental variance embodied in Eq. [5] 
However, in these experiments, —[A/p] l g is a fixed con- 
stant that need not be actually known. They therefore 
do not test the accuracy of any theoretical predictions 
for —[A/pY , i.e. the absolute degree of fractionation at 



5 



45- 
40- 
35- 
30- 
' 25- 
20- 

g 

15 
10 
5 




0.30 



0.35 
4 



0.40 

7:8 n H 
3 R 



0.45 



0.50 



FIG. 3. Predicted fractionation strength as a function of poly- 
mer reservoir volume fraction in a CP mixture with polydis- 
perse colloids, within Evans's perturbative theory applied to 
the FVAO theory, with two choices of the HS compressibility 
factor. The filled symbols show Zcs, the open ones Zwkv3 
for q = S/a = 0.4(o, .), 0.15(0, ♦), 0.1(D, ■). 



a given parental variance. Later, pre dictions made using 
the moment free energy method |45j showed reasonable 
quantitative agreement with the aforementioned experi- 
ments in terms of both the absolute degree of fraction- 
ation as well as its scaling with the parent distribution 
width; however, these were only performed at moderate 
size ratio, where the poor performance of CS has little 
impact. A combined experiment, theory and simulation 
study based on recently published data [4L| that is aimed 
at testing calculated values of — [A/p] 1 for a CP mixture 
at low q is under way. However, Fig. [3] already serves to 
illustrate that the behaviour of Zcs for short attractions 
has significant consequences not just for the overall phase 
behaviour but, perhaps more so, for supervening theories 
such as the perturbative theory of polydispersity tested 
here. 



IV. CONCLUSIONS 



showed that the known unphysical behaviour of the gas- 
liquid binodal in this regime could be attributed in part 
to the use of the Carnahan-Starling hard sphere free en- 
ergy, whose divergence at cj> = 1 allows the liquid to be 
compressed more than is physically possible at a spu- 
riously low free energy cost; although CS- independent 
failings of FVAO at low-g remain very important. Many 
other expressions for ZhSi an d details of their accuracy, 
are available We chose one of these, WKV3, for 

its similarity to CS at moderate (j) and strong difference 
at higher </> [l9j |. using it as an alternative input to the 
FVAO theory to demonstrate the sensitivity of predic- 
tions in the short attraction range regime to the choice 
of Zns, and so the importance of this choice in mak- 
ing and evaluating such predictions. In particular, we 
showed that at q < 0.2, the use of Zcs and Zwkv3 lead 
to quite different results both for the gas-liquid binodal 
and for the degree of fractionation when the system is 
polydisperse. 

In a wide variety of theories (in which CS is the de facto 
choice for Zhs) H, HH l35l - [39j . predicted gas- liquid 
binodals are found to become increasingly inaccurate as 
attraction range decreases; it seems likely that CS con- 
stitutes a specific common point of failure, in addition to 
the general difficulty in deploying a HS reference in this 
regime. We emphasise that in none of these examples is 
the unphysical divergence of CS noted, nor is any other 
choice of Zhs entertained. This could be particul arly 
important in studies such as that reported in Ref. [17| . 
where the CS reference is used alongside a Yukawa po- 
tential, a crystal phase, and a mode coupling theory. In 
such situations, a significant quantitative change in po- 
sition of the gas-liquid binodal could then qualitatively 
alter its relationship to those other elements of the phase 
diagram. Similarly, it could be important in interpreting 
findings (e.g. [36|]) in which multiple independent theo- 
ries disagreeing with simulation are in fact all dependent 
on the same CS free energy. 

Predictions of phase behaviour and related phenom- 
ena (e.g. fractionation) in systems with short-range at- 
tractions can easily be tested for dependence on Zcs by 
a similar procedure to that performed here. Otherwise, 
detailed analysis of a theory's performance may be ill- 
founded, given that large discrepancies caused by using 
the CS expression (separate to those intrinsic to the the- 
ory incorporating it) may go unnoticed and unquantificd. 



In this paper, we have studied the FVAO theory of 
CP mixtures as a representative HS reference-based the- 
ory, in the regime where the polymer to colloid size ratio, 
q = S/a m r g /a, is significantly less than unity, so that 
depletion attraction induced by the polymers between 
the colloids is short-ranged. By testing an alternative 
expression for the HS compressibility factor, Zwkv3, we 
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